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We use newly developed 3D kinetic MHD models of the heliosphere to predict 
heliospheric H I Lya absorption for various lines of sight. These predictions are 
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compared with actual Lya spectra from the Hubble Space Telescope, many of 



which have yielded previous detections of heliospheric absorption. We find that 
the absorption predicted by the models is weakly affected by both the magnitude 
and orientation of the assumed ISM magnetic field. Models with B = 1.25 — 



> . 

2.5 fiG and an angle of a = 15 — 45° with respect to the upwind direction of 
! the ISM flow generally provide the best fits to the data, but the sensitivity of 

the Lya absorption to many model input parameters makes it difficult to fully 
characterize the region of parameter space allowed by the data. We also use the 
models to assess the degree to which heliospheric asymmetries induced by the 
ISM field should be apparent in Lya absorption. An ISM field that is skewed 
with respect to the ISM flow vector results in substantial azimuthal asymmetries 
in both the hydrogen wall and heliosheath, but only the heliosheath asymmetries 
yield potentially detectable asymmetries in Lya absorption; and then only in 
downwind directions, where comparison with the data is complicated by few 
actual absorption detections and an insufficient model grid extent. 
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INTRODUCTION 



The interaction region between the solar wi n d and ambient ISM has been the subject 
of hydrodynamic modeling efforts (jParkerl Il96ll . 11963) since around the time of the first 
in situ observations of the so lar wind by Mariner 2 ( INeugebauer fc Snyderl Il962l ) and by 
Luna 2 (IGringauz et al.lll960l ). As shown in Figure 1, this interaction results in a large scale 
structure for the heliosphere that consists of three boundaries: the termination shock (TS), 
where the solar wind is shocked to subsonic speeds; the bow shock (BS), where the ISM flow 
is shocked to subsonic speeds; and in between the two the heliopause (HP), which separates 
the plasma flows of the fully ionized sol ar wind and p a rtially ionized I S M. Re v iews of the 



history of heliospheric model ing include iHolzerl (1198 
Baranov &: Izmodenovl (120061 ). 



Baranov! (11990f ). IZankl (Il999f ). and 



Voyager 1 recently encountered the TS at a distance of 94 A U from the Sun in roughly 
the upwind direction relative to the ISM flow (IStone et al.ll2005l ). However, the locations of 
the more distant HP and BS remain observationally uncertain. In general, there are very 
few observational constraints for the properties of the heliosphere beyond the TS. One of the 
few exceptions is heliospheric Lya absorption, which is observable in Hubble Space Telescope 
(HST) spectra of nearby stars. 

Unlike the ionized component of the ISM, the neutrals in the ISM can penetrate into 
all regions of the heliosphere. Charge exchange processes involving these neutrals create 
populations of hot H I that permeate the heliosphere, and it is these neutrals that produce 
absorption signatures in stellar Lyct lines observed by HST. For most lines of sight, the 
absorpti on is dominated by H I in the so-called " hydrogen wall" region in between the HP 
and BS (IBaranov et al.l Il99ll ; IWood et al.l l2005bl ). but in downwind directions absorption 
from helios heath neutrals, created by charge exchan ge between the TS and HP, can be 
paramount (jlzmodenov et al.lll999l ; IWood et al.l 120071 ). The heliospheric absorption is only 
detectable when the ISM absorption for the observed line of sight is not too broad to obscure 
the absorption. In upwind directions the interstellar H I column density (in cm -2 ) must be 
logiV(H I) < 18.2 to det ect heliospheric abs orption, but in downwind directions one must 
have logiV(H I) < 17.8 (IWood et al.l l2005al ). Astrospheric absorption from the wind-ISM 
interaction region surrounding the observed star can also sometimes be detected. 



Starting with iGayley et al.l (119971 ). there have been many attempts to use the Lya ab- 
sorption observations to test heliospheric models. The hydrodynamic models are generally 
quite successful in reproducing the observed amount of absorption, especially in upwind di 



rections where the hydrogen w all accounts for most of it (IGayley et al.lll997l ; llzmodenov et al. 
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tion of the hydrogen wall, and a validation of the models that predicted it even before it was 
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detected by HST. 



However, the exact amount of absorption predicted by the models is dependent on 
the par ameters that are assumed for the Local Interstellar Cloud (LIC) in which the Sun 
resides ( jLallement fc Bertinlll992l ). Some aspe cts of the ambient ISM are know n very well, 
such as the LIC flow speed and direction (e.g., IWitte 112004 ; iMobius et al.l 12004 ). but others 
are not known as precisely. Thus, ther e has been hope t hat t he Lya absorption can help 
constrain certain properties of the ISM. Ilzmodenov et al.l (120021 ). for example, experimented 
with numerous different models assuming different combinations of ISM proton and H I 
densities. The absorption predicted by the models does vary with the input parameters, but 
the absorption diagnostic seems to have only a modest sensitivity to most input parameters 
of interest, making it difficult to simply define a range of parameters that are consistent with 
the data. The dependence of the predicted abso rption on the nature of the hydrodynami c 
code used in the modeling is also a problem (jWood et al.l |2000| ; Ilzmodenov et al.l 120021 ). 
The source of this difficulty lies in the com plexity of H I velocity distributions within the 
heliosphere, which are non-Maxwellian (e.g.. Ilzmodenov et al.ll200ll ) and therefore can only 
be modeled with fully kinetic or complex multi-fluid codes. 

All of the models that have been compared with the data in the past have been 2- 
dimensional, axisymmetric models. Recently, 3-dimensional MHD models have become 
available that are capable of considering the effects of the ISM magnetic field on heliospheric 
structure, while still maintaining a sufficiently sophisticated treat ment of the neutrals to 



properly consider them and the plasma in a self-consis tent manner (Ilzmodenov et al.l 12005 



Izmodenov fc Alexashovl 120061 ; iPogorelov fc Zankll2006l ). Figure 1 presents results o f calcu - 
lations made with a 3D kinetic MHD model of the heliosphere by Ilzmodenov et al.l ( 20051 ). 
It shows that the inclusion of even a modest ISM field can indeed affect the shape of the 
global heliosphere. We will determine here whether this also has significant effects on the 
Lya absorption. The nature of the magnetic field in the ISM immediately outside the he- 
liosphere is poorly known, so we will also assess the sensitivity of the Lya absorption to 
changes in the assumed ISM field strength and orientation. In doing so, we consider many 
more HST-observed lines of sight than have been used in prior data-model comparisons. 



2. THE CHOSEN SAMPLE OF HST Lya OBSERVATIONS 

The amount of heliospheric Lyct absorption depends greatly on the direction of the 
observed line of sight. The greatest spatial dependence is on the poloidal angle 6 between 
the line of sight and the upwind direction of the ISM flow. Clearly it is advantageous to 
consider many different lines of sight with a wide variety of #'s in comparing the heliospheric 
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absorption predicted by models with the data. Considering a variety of directions is even 
more important when testing 3D MHD models, which can yield heliospheric structures that 
are not axisymmetric and therefore will have absorption predictions that are dependent on 
the azimuthal angle as well as being dependent on 6 (see Fig. 1). Past data-model com- 
parisons considered no more than six HST-observed lines of sight, which indiv idually have 



eithe r provided real detection s of heliospheric absorption or merely upper limits (IWood et al. 



2000l ; llzmodenov et al.ll2002l ). This is a rather small number of lines of sight even for testing 
axisymmetric models, let alone the 3D ones. Howe yer, the number of heliosp heric absorp- 
tion detections has recently increased significantly (IWood et al.l l2005bl . 120071 ). so it is well 
worthwhile to reassess the sample of available HST data to select a larger sample of spectra 
to test the 3D kinetic MHD models. 



Wood et al.l (l2005bl ) provide a complete list of HST-observed Lya spectra that are ap- 
propriate for our purposes, all of which have been analyzed to measure ISM H I column 
densities, to search for evidence of heliospheric/astrospheric absorption, and to measure 
stellar Lya fluxes corrected for the contaminating ISM absorption. Figure 2 is a sky map 
in ecliptic coordinates of lines of sight with a stellar Lya line observed by HST. All these 
spectra have sufficient spectral resolution to permit a reasonably precise search for helio- 
spheric absorption. The boxes indicate the 11 lines of sight that actually yield detections of 
heliospheric absorption. All the other lines of sight yield nondetections. 

Many of the detections are clustered around the upwind direction of the ISM flow. 
The advantageous nature of upwind lines of sight for detecting heliospheric absorption is 
consistent with model predictions, which suggest that the deceleration of H I in the hydrogen 
wall relative to the ISM flow should be largest in these directions. This results in a greater 
separation of the heliospheric absorption from that of the ISM, thereby m aking it easier to 
detect heliospheric absorption in upwind lines of sight (IWood et al.ll2005bl ). 



There is also a cluster of three detections very c lose to the downwin d direction. Initial 
analysis of these Lya spectra did not yield detections (IWood et al.ll2005bl ) . However, we have 
found that the stellar Lya profiles reconstructed for 9 > 160° lines of sight are systemati- 
cally blueshifted from the stellar rest frames, in dicating the presen ce of very broad, shallow 
absorption on the red side of the Lya profiles (IWood et al.l 120071 ). This is the exactly the 



sort of absorption signature one expects from heliosheath neutrals (as opposed to hydrogen 
wall neutrals). Since very downwind lines of sight looking down the tail of the heliosphere 
will have very long path lengths through the heliosheath, it is in the most downwind lines 
of sight where one might expect to see this broad absorption. Thus, we now consider these 
three lines of sight to have detections of heliospheric absorption, though the nature of these 
detections is rather different from the others. 
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Our goal is to select a sample of HST-observed lines of sight from Figure 2 to use 
for comparing observed and predicted heliospheric Lya absorption. Obviously we start 
by choosing the 11 detections, which actually provide quantitative measurements of the 
absorption. We add to these detections nine nondetections (diamonds in Fig. 2) that at 
least provide upper limits for the amount of absorption that might be present in those 
directions. These nondetections are chosen to sample parts of the sky not covered by the 
detections. Another major selection criterion is ISM H I column density. Lines of sight with 
low ISM column densities are preferable since they provide more restrictive upper limits on 
heliospheric absorption. Data quality (i.e., resolution, signal-to-noise) also plays a role in 
choosing which nondetections to consider. 

Numbered symbols in Figure 2 indicate the final sample of 20 lines of sight to be used in 
our data-model comparisons. The stellar identifications of these lines of sight are indicated 
in Figures 3 and 4, along with the Lya spectra, which are displayed in order of increasing 
9. We focus only on the red side of the Lya absorption profile, where the heliospheric ab- 
sorption resides. The fluxes are normalized to the intrinsic st ellar Lya profile reco nstructed 



in the original analysis of the data. We refer the reader to IWood et al.l (]2005bl ) and ref- 
erences therein to see the full Lya spectra and descriptions of their analysis. The dotted 
green lines in the figure show only the ISM absorption based on these analyses. For the 
heliospheric absorption detections, there is excess absorption observed beyond that from the 
ISM. Successful heliospheric models should predict the right amount of excess absorption to 
fit the data for these lines of sight. For the nondetections, the ISM absorption fits the data 
reasonably well. In these cases, successful heliospheric models should predict essentially no 
significant absorption beyond that from the ISM. 

The three 9 > 160° detections (#18-#20 in Figs. 2-4) are special mentioned 
above. The original reconstructed stellar Lya profiles suggest no heliospheric absorption, but 
the blueshifts of these profiles away from their stellar rest frames implies that these profiles 
are inaccurate. We can infer the amount of heliospheric absorption in these directions by 
constructing a stellar profile forced to be centered on the stellar radial velocity and then 
seeing how much of the red wing of that profile must be absorbed to yield the original profile. 
The shaded regions in Figure 3 for these three downwind lines of sight indicate this excess 
absorption, where the uncertainties are estimated by allowing the stellar radial velocity to 
be ±3 km s _1 from its measured value. This excess absorption cannot be extended to lower 
velocities closer to the center of the Lya line because near line center stellar Lya profiles 
cannot be assumed to be symmetric and centered on the stellar rest frame. Stellar Lya 
profiles often have s elf-reversals near li ne center, which are often asymmetric. For details 



about all of this, see IWood et al.l (120071 ). The important thing to note here is that for these 



three very downwind lines of sight, the absorption predicted by the models should not fit 
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the data but should instead fall within the shaded regions. 

Finally, the requirement that the intrinsic stellar Lyct profile be within ±3 km s _1 of 
the stellar rest frame allows us to compute upper limits for the amount of broad heliosheath 
absorption that can be present for lines of sight without detected heliospheric absorption. 
Thick dashed lines in Figure 3 show these upper absorption limits, but only for downwind 
lines of sight (9 > 110°) where the broad heliosheath absorption is potentially prominent. 
Absorption predictions from the models must lie above these limits to be consistent with the 
data. The dashed lines cannot be extended to low velocities close to line center for the same 
reason that the shaded regions of the 9 > 160° lines of sight are limited to the wings of the 
line (see above). 



3. THE INTERSTELLAR MAGNETIC FIELD'S EFFECTS ON Lyct 

ABSORPTION 

Figures 3 and 4 compare the HST Lya data with the helio spheric absorption pre- 



dicted by 3D kinetic MHD m odels of the heliospheric interface ( llzmodenov et al.l 12005 



Izmodenov fc Alexashovl 120061 ). assuming various directions an d magnitudes for th e ISM 



field. The models used here are of the type initially developed by lBaranov fc Malamal (11993 



19951 ). with a fully kinetic treatment of neutral hydrogen within the heliosphere to pro- 



vide the most precise c omputations of the velocity distribution functions of the neutrals. 



Izmodenov et al.l (120051 ) expanded the 2D axisymmetric Baranov & Malama code to a fully 
3D geometry, and also added the capability of including an interstellar magnetic field in 
the model. The code separates all heliospheric H atoms into several populations: 1. origi- 
nal interstellar atoms and other atoms originating outside of the bow shock, 2. secondary 
interstellar atoms originating between the bow shock and heliopause, 3. atoms originating 
between the heliopause and termination shock, and 4. atoms originating in the supersonic 
solar wind. We calculate number densities, temperatures, and bulk velocities for these pop- 
ulations along the lines of sight toward the observed stars by taking moments of the velocity 
distributions. The heliospheric absorption for each line of sight is computed from these 
traces of density, temperature, and flow velocity. With this methodology, we are making 
the approximation that the velocity distribution functions of the individual populations are 
locally Maxwellian. 
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3.1. Previous Constraints on the Local ISM Field 



The nature of the interstellar magnetic field surrounding the Sun is poorly known, 
though some observational constraints exist. The global Galactic field has a magnitude 
of 1.6 ± 0.2 /xG and is direct ed towards a Galactic longitude of I = 96 ± 4°, but there is 
substantial local variability (IRand &: Kulkarnil Il989l ). meaning that the actual local field 
could be significantly higher or lower and could be in a completely different direction. 

A ~ 4° discrepancy exists between the flow vectors of interstellar He and H within 
the solar system, and the most promising explanation for this is that the LIC's magnetic 
field is skewed with respect to the ISM flow seen by the Sun, which can deflect the flow 
of interstellar hydrogen atoms in the heliosphere (jLallement et al. 2005 ). Helium atoms 
are not affected in this manner, since their charge exchange cross sections are much lower 
than hydrogen and they are therefore effectively blind to the presence of the heliosphere. If 
correct, this interpretation identifies a plane in which the ISM field mus t lie, which h a ppens 



to be inconsistent with the orientation of the global Galactic field. lOpher et al.l (120061 ) 
have argued that an ISM field that is a = 30° — 60° from the apparent flow direction can 
potentially explain Voyager 1 and 2 observations of energetic particles flowing inwards from 
the TS. The satell ites have both seen these particles, but flowing in opposite directions. 
Opher et al.l (120061 ) demonstrate that asymmetries in the heliospheric structure induced by 
a skewed ISM magnetic field can potentially cause this effect, thanks in part to the satellites' 
positions on opposite sides of the ecliptic plane. 

A magnetic field much stronger than the global Galactic field has been proposed to 
explain an apparent pressure imbalance between the hot, ionized plasma that dominates the 
Local Bubble and the warm, partial ly neutral clouds that lie within it. The Local Bubble, 
the cavity in which the Sun resides ( Lallement et al. 2003h. is belie ved to account for much 
of the soft X-ray background radiation (see also ISmith et al.l 120061). The se X-rays seem to 
suggest thermal pressures of P/k ~ 15,000 cm -3 K ( jSnowden et al.lll998l ). In contrast, the 
LIC and other warm clou ds within the Local Bubble appear to h ave much lower pressures 
of P/k ~ 2280 cm" 3 K Jjenkinsl liooi iRedfield fc Linskvi 120041 ). Such a large pressure 
imbalance within the local ISM seems unlikely. 

One way out of this dil emma is to propose t hat the LIC is supp orted by a strong 
magnetic field of order 7 /xG ( Cox fc Helenius 2003 ; Florinski et al. 20041 ). This would seem 
to be disallowed by heliospheric models, which imply that such large magnetic pressures 
would force the termination shock well i nside the 94 AU distance measured by Voyager 1 
( iGloeckler et al.l 119971 ; IStone et al.l 120051 ) . The only way to avoid this problem is for the 
field to be nearly parallel to the ISM flow. However, this would contradict the evidence 
mentioned above that the LIC field is not parallel to to the LIC flow (i.e., the discrepant H 
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and He flow vectors within the solar system, and the opposite flow directions of termination 
shock particles observed by Voyager 1 and 2). Thus, large fields are still difficult to reconcile 
with heliospheric observations and models. An alternative solution to the Local Bubble 
pressure problem is that the pressures normally estimated from the soft X-ray emi ssion are 
too h i gh for various reasons: contamination from heliospheric foreground emission (j Cravens 



2000l ; iLallementl l2004al ; iKoutroumpa et al. 



the walls of the Local Bubble (Lallement 



2006), contamination from X-ray emission from 



2004b|), and improper assumption of collisional 



ionization equilibrium (IBreitschwerdtl 1200 ll ) . 



3.2. Absorption Constraints on the ISM Field 

Most of our models are computed assuming a modest magnetic field of B = 2.5 /iG, 
which can compress the TS somewhat depending on the field orientation, but not enough to 
be inconsistent with the Voyager 1 encounter distance of 94 AU (see Fig. 1). We as sume that 



the magnetic field is oriented within the plane suggested by lLallement et al.l (120051 ) . Figure 3 
shows the absorption predicted by models with different field orientations within this plane, 
for angles ranging from parallel to the ISM flow (a = 0°) to perpendicular to the ISM flow 
(a = 90°). Table 1 lists the ecliptic and Galactic coordinates that correspond to these field 
directions. Figure 4 shows the absorption predictions for three models with a = 45°, but 
with different field strengths of 0, 1.25, and 2.5 /iG. The ISM hydrogen and proton densities 
assumed in these models are n^H I) = 0.18 cm -3 and n 0O (H + ) = 0.06 cm~ 3 , the ISM flow 



speed is Voo = 26.4 km s x , and the temperature is Too = 6400 K (jlzmodenov et al.ll2005l ). 



Figure 3 shows that in upwind directions, where most of the heliospheric detections lie, 
the absorption decreases with increasing a. To better illustrate this behavior, Figure 5 plots 
the absorbed Lya flux predicted by the Figure 3 models versus a. No one model clearly 
fits the data better than all the others, presumably due to the absorption dependence on a 
being only a modest one. This also may be indicative of the systematic uncertainties in the 
estimation of the ISM absorption. The 61 Vir line of sight seems particularly discrepant. 
This is probably due to difficulties in remo ving geo coronal em ission blended with the red 



side of the Lya absorption line (see Fig. 1 in lWood et al.ll2005bl ). so this line of sight should 
be regarded with caution. 

The a = 60° and a = 90° models underestimate the absorption in all upwind directions 
(see Figs. 3 and 5), so perhaps these models should be considered inconsistent with the data. 
If one ignores the problematic 61 Vir line of sight, the a = 0° model is a worse fit to the 
data than the a = 15° — 45° for all but the 36 Oph line of sight. Considering the evidence 
described above that the ISM field is skewed from the flow direction, a = 0° seems unlikely 
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anyway. Thus, we consider a = 15° — 45° to represent the most likely field or i entati ons for 



the local ISM, which overlaps the a = 30° — 60° range quoted by lOpher et al.l (120061 ). 



These conclusions are based on B = 2.5 fiG models, but Figure 4 shows that assuming 
B = 1.25 fiG does not change the absorption very much. Thus, at least for these low- 
to-moderate field strengths our conclusions concerning a are relatively sound. However, 
Figure 4 shows that decreasing the field all the way to zero does change the absorption, 
with the absorption being somewhat higher in upwind directions. The greater amount of 
absorption upwind is somewhat surprising considering that the hydrogen wall is narrower 
for B = (see Fig. 1). However, the inclusion of even a modest ISM field weakens the bow 
shock and lowers the H I density in the hydrogen wall, more than offsetting the broader 
width of the wall. 

It should be stated that definitive conclusions are difficult to make at this point since the 
Lya absorption is at least somewhat sensitive to other input parameters, su ch as the assumed 



ISM densities and temperature that also have observational uncertainties (jlzmodenov et al. 



20021 ). though not nearly as large as those involving the magnetic field. A time-consuming 
thorough exploration of parameter space would be necessary to fully characterize the con- 
straints on ISM properties provided by the Lya data. 

We have focused so far on comparing the models with data in upwind directions, where 
absorption from the hydrogen wall is dominant. Any conclusions drawn from downwind 
directions, where heliosheath absorption is more prominent, will be far more tentative. The 
primary reason for this is that the grid used for our current models only extends 500 AU 
from the Sun, which is not nearly far enough to capture all of the heliosheath absorption for 
6 > 120°. Thus, the absorption predictions shown in Figures 3 and 4 for these directions 
will underestimate the amount of absorption that the models would really predict if the grid 
were extended further downwind. 

Many of the models seem to predict too much absorption downwind even with the lim- 
ited grid extent, particularly in the velocity range of 80 — 120 km s" 1 . But there is another 
potential difficulty with downwind absorption that concerns the treatment of the plasma 
in the models. Although a fully kinetic treatment is applied to the neutrals, the plasma 
is assumed to be a single Maxwellian fluid throughout the heliosphere. However, this is 
clearly a poor approximation, since pickup ions, for example, have non-Maxwellian velocity 
distributions and are no t thermalized wi t h the solar wind inside the termination shock (e.g., 



Gloeckler fc Geisd 12004 ) . iMalama et al.l (120061 ) have replaced the simple single- fluid plasma 



treatment in the 2D Baranov & Malama code with a complex multi-component representa- 
tion of the plasma. This more sophisticated plasma treatment does not result in significantly 
different hydrogen wall absorption in upwind directions, but we have found that is does result 
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in a significant reduction in heliosheath absorption which can potentiall y alleviate the prob - 
lems these models have in predicting too much downwind absorption (jWood et al.l 120071 ). 
The most meaningful comparison with the data in downwind directions would therefore re- 
quire that our 3D models also utilize such a multi-component plasma treatment, as well as 
having grids that extend far enough downwind to capture all the heliosheath absorption. We 
leave such computationally intensive modeling for a future paper. 

We note that asymmetries in the heliospheric structure induced by the ISM field are 
evident in the downwind absorption predicted by the models. For example, there is significant 
model dependence in the absorption towards DK UMa (9 = 112°), but no significant model 
dependence towards r Cet (9 = 123°). Since 9 is similar for these stars but their azimuthal 
angles quite different (see Fig. 2), this difference in behavior must be due to azimuthal 
variability, which can only be due to magnetic field induced asymmetries (see §3.3). 



Finally, we note that iPogorelov fc Zankl (120061 ) have also developed 3D MHD helio- 
spheric models that treat neutrals in a self-consistent manner with the plasma. This code 
uses a less sophisticated 2-fluid treatment for the neutral H velocity distributions, but it 
includes the effects of the interplanetary magnetic field on heliospheric structure as well 
as the ISM field. Absorption has been computed for a limited number of these models, 
yielding results qualitatively similar to th ose reported here, with a bsorption decreasi n g wit h 
both a and the magneti c field strength (jWood et al.l 120061 ) . Like llzmodenov et al.l (120051 ) , 
Pogorelov Zankl (120061 ) demonstrate that the 3D MHD models ca n potentially reproduc e 
the shift between the H and He flows observed within the solar system (ILallement et al.ll2005l ). 
though they emphasize that the magnitude of the shift depends not only on the strength 
and orientation of the magnetic field, but also on the ISM neutral hydrogen density. 



3.3. Quantifying Expected Absorption Asymmetries 



Figure 1 illustrates that even a modest ISM field can result in a heliospheric structure 
that is signi ficantly asymmetric, consistent with other mode ls that also predict asymmetries 
of this sort (IRatkiewicz et al.lll998l ; IPogorelov fc Zan kl l2006h . In the bottom half of Figure 1, 
the heliosheath (between the TS and HP) is narrower than in the upper half, but the hydrogen 
wall (between the HP and BS) is wider. One might imagine that this would result in 
corresponding Lya absorption asymmetries. In other words, the Lya absorption should be 
azimuthally dependent as well as 9 dependent. 

Some evidence that the models do indeed predict azimuthally dependent absorption 
in downwind directions is mentioned in §3.2. However, in order to properly quantify the 
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degree of absorption asymmetry expected based on the models, it is necessary to compare 
absorption predictions for lines of sight with identical values but different azimuthal angles 
(0). This cannot be easily done with the set of observed directions in Figures 3-4, which are 
scattered randomly about the sky 

Thus, in Figure 6 we show the heliospheric Lya absorption predicted by the B = 2.5 /xG, 
a = 45° model for various angles, with 9 fixed in the five panels of the figure. The azimuthal 
angle is defined such that = 0° and = 180° are in the plane of the ISM magnetic field. 
This is the plane in which the heliospheric structure is portrayed in Figure 1. The = 0° 
direction would be associated with the upper half of Figure 1, and = 180° would be 
associated with the lower half. 

In general, the hydrogen wall will be responsible for the steep, saturated absorption edges 
of the absorption profiles in Figure 6, which are particularly prominent in upwind directions 
(e.g., located at 20-30 km s _1 in the 9 = 30° panel), while the heliosheath is responsible for 
the broad, unsaturated absorption wings that extend to high velocities, which become more 
prominent in downwind directions. Very little dependence is apparent in upwind directions. 
This is a bit surprising given the hydrogen wall asymmetries apparent in Figure 1. However, 
it turns out that azimuthal density variations in the hydrogen wall offset the azimuthal width 
dependence. For example, although the hydrogen wall is narrower for = 0° (corresponding 
to the upper half of Fig. 1) than for = 180° (corresponding to the lower half of Fig. 1), this 
is offset by higher hydrogen wall densities in the = 0° direction, so integrated H I column 
densities are actually not very different. 

Figure 6 shows that for 9 = 60 — 90°, the particularly narrow hydrogen wall at = 0° 
results in somewhat less hydrogen wall absorption than other directions, but the difference is 
so small it would be very difficult to detect in practice. In contrast, the heliosheath is thicker 
at = 0° and Figure 6 shows that this does in fact lead to the broad heliosheath absorp- 
tion being significantly stronger in this direction than others, though it is only in downwind 
directions (9 > 90°) where this azimuthal dependence becomes potentially detectable. As 
discussed in §3.2, actual comparisons with the data are currently problematic in these direc- 
tions. The asymmetries seen for the 9 = 90° and 9 = 120° surely extend to 9 = 150° as well, 
but Figure 6 does not show this due to the limited grid extent of the models. According to 
Figure 6, the heliosheath absorption is at a minimum in directions normal to the plane of 
the ISM field (0 = 90° and = 270°), indicative of magnetic compression of the heliosphere 
in those directions, leading to shorter distances through the heliosheath for those lines of 
sight. 



4. SUMMARY 



We have compared H I Lya absorption profiles predicted by 3D kinetic MHD models 
of the heliosphere with a large selection of Lya lines observed by HST, including many lines 
of sight with detected heliospheric absorption. The primary purpose of this comparison is 
to assess the sensitivity of the predicted absorption to changes in the ISM magnetic field 
properties assumed in the model. Our results are as follows: 



We find that the H I Lya absorption has a modest sensitivity to both the strength and 
orientation of the ISM magnetic field. Focusing on upwind directions where most of 
the HST detections of heliospheric absorption reside, the models presented here with 
B = 1.25 — 2.5 fiG and a = 15 — 45° appea r to provide the best fits to the data, 



consistent with constraints from other sources (Glocckler et al. 1997; Izmodenov et al. 



20051 : lLallement et al J 1200.4 lOpher et alJl200fih . 



2. However, since the Lya absorption is sensitive to other model input parameters, such 

as the ISM H I density, which have not been varied here, the region of parameter 
space that yields acceptable fits to the data will be complex. It will be very difficult, 
perhaps impossible, for the Lyct absorption by itself to yield a unique set of model input 
parameters that fit the data. Nevertheless, the dependence of the absorption on many 
ISM parameters means that the absorption does provide one constraint on heliospheric 
models that is worthy of consideration in assessing how precisely the models reproduce 
reality. 

3. The models show that an ISM field that is skewed with respect to the ISM flow vector 

results in substantial azimuthal asymmetries in the heliospheric hydrogen wall. Sur- 
prisingly, these asymmetries do not result in significant asymmetries in Lya absorption 
from the hydrogen wall, since densities within the wall vary in such a way as to cancel 
out the effects of the spatial asymmetries on hydrogen wall column densities. 

4. The models also show that a skewed ISM field results in significant azimuthal asym- 

metries in the heliosheath, and unlike the hydrogen wall these asymmetries do yield 
significant azimuthal absorption dependence, at least in downwind directions where the 
heliosheath absorption is prominent. These directions are clearly the best places to look 
for azimuthal dependences in Lya absorption, but there are problems with doing this 
in practice. One is simply that we have few downwind detections of heliospheric Lyct 
absorption. Another is that the heliosheath absorption that dominates in downwind 
directions should ideally be modeled using a complex multi-fluid plasma treatment. 
And finally, the model grid must be extended a much longer distance downwind than 
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present models to capture all the heliosheath absorption in these directions. We hope 
to perform such computationally intensive modeling in the future. 
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Table 1. Model ISM Field Orientations 



a (deg) 


Ecliptic 


Coord. 


Galactic Coord. 




l e (deg) 


b e (deg) 


1 (deg) 


b (deg) 
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Fig. 1. — The locations of the termination shock (TS), heliopause (HP), and bow shock 
(BS) for a model including a B = 2.5 /iG ISM field (solid lines), and a model with no ISM 
magnetic field (dashed lines). The directions of the LIC flow vector (Vlic) and magnetic 
field (Buc) are indicated. The distance scale is in AU. The region between the HP and 
BS is sometimes called the "hydrogen wall" and the region between the TS and HP is the 
"heliosheath." 
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Fig. 2. — Sky map in ecliptic coordinates of all HST-observed lines of sight with useful Lya 
spectra. The numbered symbols indicate spectra that we will compare with model predictions 
of Lya absorption (see Figs. 3 and 4). Boxes indicate lines of sight with detected heliospheric 
absorption. The plus signs and diamonds are both lines of sight with nondetections of 
heliospheric absorption. The diamonds indicate lines of sight selected to provide upper 
limits for absorption in those directions. The filled and open circles indicate the upwind and 
downwind directions of the local ISM flow vector, respectively. 
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Fig. 3. — The red side of the H I Lya absorption line (histogram) for the selected stars from 
Fig. 2, where the stars are placed in order of increasing angle from the upwind direction 
of the ISM flow (#). In each panel, the dotted green line is the ISM absorption alone. 
Absorption predictions are shown for heliospheric models computed assuming six different 
ISM field orientations, as quantified by a, the angle between the field and the ISM flow 
direction (see the 36 Oph panel for line identifications). For many downwind lines of sight 
(9 > 110°), dashed lines show upper limits to the amount of absorption that can be present 
- absorption predictions from the models must lie above these limits to be consistent with 
the data. For the three most downwind lines of sight, the shaded regions indicate the amount 
of absorption that the models should predict if the real stellar Lya profile is centered on the 
stellar rest frame rather than blueshifted as suggested by the original recontructed profile. 
For these lines of sight, the absorption predicted by the models should not fit the data (which 
are dotted histograms in these cases) but should instead fall within the shaded regions (see 
§2). 
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Fig. 4. — A figure analogous to Fig. 3, but in this case the absorption predictions are for 
three a = 45° models that assume different ISM magnetic field strengths (see the 36 Oph 
panel for line identifications). 
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Fig. 5. — For six upwind lines of sight, the predicted wavelength-integrated Lya flux absorbed 
by heliospheric H I beyond that absorbed by the ISM is computed for the six models from 
Fig. 3 and plotted versus a, the ISM field orientation relative to the ISM flow direction. The 
fluxes are normalized to the observed amount of flux absorbed, so in each panel a flux of 1 
(dashed lines) corresponds to agreement with the data. 
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Fig. 6. — An illustration of the directional dependence of H I Lya absorption predicted by a 
3D MHD heliospheric model assuming an ISM field strength and orientation of B = 2.5 /xG 
and a = 45°, respectively. Absorption is shown for five values of the poloidal angle 9 (the 
angle between the line of sight and the upwind direction of the ISM flow), and eight values 
of the azimuthal angle (where the plane of the ISM magnetic field is in the = 0° and 
= 180° directions). The model grid does not extend far enough downwind to properly 
search for azimuthal absorption variations at 9 = 150°. 



